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ABSTRACT 


The power spectral techniques of generalized harmonic analysis 
were applied to determine the phase and amplitude of the transfer 
function between the vertical component of turbulence and the 
resultant pitch rate of a small swept wing jet fighter. The power 
spectrum analysis was based on actual flight test data. In addition 
the amplitude and phase of the aircraft's transfer function are 
determined theoretically by:.two different methods and compared with 
the results of the power spectrum analysis. It was established that 
the power spectrum method, as applied to random turbulence, can be 
used to determine an aircraft's response to this random input on a 
continuous basis and that power spectrum techniques have considerable 


potential value in the study of aircraft gust loading. 
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INTRODUCTION 


Atmospheric turbulence and the resulting aircraft disturbances 
from stabilized flight due to this turbulence are random and thus 
aperiodic and non-dissipative. These disturbances can only be de- 
fined on a non-precise basis with certain probability distributions 
and average values since they are random in nature. This leads to 
the application of probability and statistical concepts, which, in 
combination with classical Fourier techniques, form a general 
approach for the analysis of problems of this type. This technique 
is termed the power spectrum method. 

Power spectral methods were derived from economic time studies 
at the beginning of this century. The initial applications were also 
directed toward turbulence studies and the Brownian motion of parti- 
cles. During World War II power spectral methods were used to analyze 
problems in servomechanisms and communications; however, the applica- 
tion of power spectral analysis to the field of aeronautics, and in 
particular to the study of the dynamic response of aircraft to tur- 
bulence, is relatively new and still in the early stages of development. 

Other recent studies in the field of aeronautics involving 
power spectral techniques have been applied to: 

a. Wind tunnel turbulence 

b. Gust loads on aircraft 

c. Acoustic Noise 

d. Buffeting loads 

e. Electronic noise in guidance atid gunnery systems 

In particular, considering the gust loading and dynamic response 
of aircraft, statistical inference is a realistic approach since the 
problem is essentially one of a random process phenomena, Perhaps a 
more meaningful presentation of the effectiveness of the power spec- 
trum method can be made if the earlier discrete gust design methods are 


mentioned in comparison. 
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The earliest discrete gust approach involved the assumptions of 
a rigid airplane, a ramp type gust, and a specified maximum vertical 
velocity. Subsequent discrete gust design methods were adjusted for 
different wing chords and pitch characteristics between different air- 
craft by expressing data in terms of the mass ratio, 2 - The dis- 
crete gust concept was based primarily on the premise that the rela- 
tive loads for single isolated gusts are a measure of the relative 
loads in a sequence of gusts. 

However, in recent years a more general approach to turbulence 
analysis has been needed in order to account for the greater variety’ 
and degree of aircraft flexibility, mass distribution, configuration 
and dynamic response characteristics, as well as the varied operating 
speeds and altitudes. Therefore, the power spectral techniques of 
generalized harmonic analysis were developed for application to tur- 
bulence and dynamic response studies. 

As conditions changed and the variety of aircraft and operating 
altitudes/airspeeds increased, the concept of the discrete gust design 
method was broadened to include flexible bodies as well as rigid 
bodies. Although the flexible aircraft was mathematically flown 
through discrete gusts of varying wavelengths, the discrete gust con- 
cept did not adequately account for the continuous nature of atmose- 
pheric turbulence, particularly as applied to lightly damped dynamic 
systems. Latest design philosophy embodies both the discrete gust 
concept and the power spectral approach, with the design initially 
analyzed and checked using the discrete gust method and then examined 
in greater detail using the power spectrum analysis. Further develop- 
ment of the power spectral method may result in its use in the future 
entirely for atmospheric turbulence studies. 

The main advantages of the power spectral approach to turbu- 
lence analysis are outlined in Ref. lL: 


1. It provides a more realistic representation of the «=~ - 
continuous nature of atmospheric turbulence. 

2. It allows airplane configurations and response character- 
istics to be taken into account in a rational manner. 

3. It allows more rational consideration of design and 
operational variations such as configuration changes, mission 
changes and airplane degrees of freedom. 
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Most of the present day applications of power spectral methods 
are concerned with the analysis of the magnitude and distribution of 
the random data; for example, the intensity of atmospheric turbulence 
and its variation with frequency. However this thesis considers a 
more unique application of power spectral techniques by establishing 
the transfer function between pitch rate and turbulence for the 
F-100C airplane based upon the analysis of a random input, sho the 
vertical component of atmospheric turbulence. This application goes 
one step beyond the analysis of random data by considering its effect 
upon the dynamic response of an airframe, The particular transfer 
function, = , 1s considered to provide a basic indication of the 
F-100 airplane's dynamic response characteristics in a turbulent en- 
vironment. To measure the accuracy of the results obtained through 
the power spectral analysis, a comparison between 2. » 48 deter=- 
mined from this method, is made with the values of a » as deter- 
mined by solving the longitudinal equations of dynamic motion for 
the F-100C airplane. 
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DISCUSSION OF THE CHARACTERISTICS AND ANALYSIS OF RANDOM DATA 


Any recorded data which describe a physical process can be 
broadly classified into either one of two categories: deterministic or 
non-deterministic. Deterministic data are information which describe a 
process which can be represented and predicted with reasonable accuracy 
by exact mathematical relationships. Some common examples of this type 
of process are the electrical potential across a condenser as it dis- 
charges through a resistor or the motion of a satellite in a known 
orbit around the earth. Examples of non-deterministic data are the 
variation of an aircraft's vertical velocity in turbulence and the 
variation in the magnitude of the lift developed on an airfoil encounter- 
ing buffeting in rough air. Due to the inherent randomness of the 
turbulence this non-deterministic data cannot be described, nor can its 
behavior be predicted by explicit mathematical relationships. In 
practical terms, the decision as to whether or not physical data are 
classified as being deterministic or non-deterministic (random) is 
based upon the ability to reproduce the data repeatedly through con- 
trolled experiments. If the experiment or phenomena which produced 
the data can be repeated with the same results, within the tolerances 
of experimental error, then the data obtained are classified as deter- 
ministic. Conversely, if identical results cannot be obtained upon 
repeating an experiment, then the data are considered to be non- 
deterministic or random in nature. 

In considering random data, therefore, it is expected that an 
particular observation producing this data will effectively represent 
only one of many possible outcomes which might have occured as a result 
of this experiment. For example, assume that the lift developed on an 
airfoil during buffeting in flight is recorded as a function of time. A 
specific lift time history will be obtained as shown in Fig. 1. 
However, if a second airfoil of identical shape is flown through a 
similar type of disturbance simultaneously, it is likely that a different 
lift time history would result. Therefore, the lift time history for 
any one airfoil encountering a particular finite disturbance is only 
one example of an infinitely large number of time histories which might 


have occurred. If the random buffeting is observed over a finite time 


18 





LIFT 


TIME 
te 
- 
¢ a 
€ 
L4FET 
TIINE 
v4 
d,, 
LIFT 
TILNE 
M 
Lo 


Figure 1. Sample records showing the lift developed on an airfoil 
in turbulence, 
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interval and is graphically or mathematically represented by a single 
time history, then this time history is referred to as a sample record 
by Bendat and Piersol (2). The collection of all possible sample 
records of lift versus time which the random buffeting might have pro- 
duced is called a random process or stochastic process in Ref. 2. 
As a result, a sample record of observed data for a random physical 
occurrence can be considered to be only one physical outcome of a 
random process. 

The properties of a physical experiment or occurrence which 
are described by a random process can be determined from the average 
values of the experimental data as computed over the range of sample 
records obtained. A collection of géample functions which comprise a 
random process is illustrated in Fig. 2. By taking the instantaneous 
values of each of the sample functions at time ti» summing the values 
and dividing by the number of sample functions, the mean value of the 


random process for a particular time t, may be computed. In addition, 


1 
a joint moment or correlation between the values of the random pro- 

cess at two different times can be computed by taking the average of 

the product of the instantaneous values at two different times ti) and 

th + J. The product obtained in the manner described above is termed 
the auto-correlation function. The mean value, ~/, a.) , and the auto- 


correlation function (t,,¢,77) are defined as: 


A 
Odie) Lam 2, y_&) (1) 


Al 
/ 
dQ (t,t, 4+Z)=hw HD, X(t) +X, (t, +7) (2) 
V-> 90 4=/ 
According to Ref. 2, the random process is stationary if 
YK (t,) and*’y (t,, is +7) do not vary as time t is varied. Ina 


more general sense, when the properties of a random process, as com- 
puted over short intervals, do not vary significantly from one interval 


to the next, then the process is considered to be stationary. 
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Figure 2. Ensemble of sample functions forming a random process. 
(Reproduced from Ref. 2). 
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Applying the concept of stationarity to atmospheric turbulence, 
it has been shown in Ref. 1 that the conditions of stationarity and 


homogeneity apply to large scale patterns of air motion within regions 


of one hundred miles and time durations of one hour. In cases of severe 


turbulence, such as encountered in a thunderstorm, these time and dis- 
tance estimates would not be valid. However, the time duration for 
homogeneous and stationary conditions may still be almost five minutes 
in length. It is usually possible to determine the properties of a 
stationary random process by computing time averages over specified 
sample functions in the ensembles (collection of sample functions). 
Consider the pth sample for the random process illustrated in 


Fig. 2. The mean value, Ahr (A), and the autocorrelation function, 


JE Ca 4) , for the yeh sample function, according to Ref. 2,are: 


An (&) = him. $/ Ke (Oot — @) 


~ 
Re(qh)=Lim $f Kg (OiXglrDdt 
Random processes may be classified as shown in Fig. 3, repro- 


duced from Ref. 2, as being either stationary or non-stationary, and 
stationary random processes may be further subdivided as either 
ergodic or non-ergodic. The random process is ergodte if it is sta- 
tionary and Ly (4) and, (7 )do not vary when computed over 
different sample functions. Therefore, all properties of ergodic ran- 
dom processes can be determined by computing over a single sample 
function. In actual physical phenomena, as outlined in Ref. 2, most 
random data, including atmospheric turbulence, are ergodic so that 
the properties of the data can be computed from time averages of 
individual sample records. 

Two types of statistical functions are used to describe the 
basic properties of random data: (a) the autocorrelation function 
and (b) the power spectral density function. The autocorrelation 


function and the power spectral density function provide an analysis 
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Figure 3. Classifications of random data (Reproduced from Ref. 2.) 
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of the amplitude of the random data with respect to time and frequency 
respectively. These functions are Fourier transform pairs, and there- 
fore provide an analysis of the data in two different formats. | 

The autocorrelation function for random data provides an indi- 
cation of the general dependence of the amplitude of the data at a 
particular time on the value of the data at a different time. An 
autocorrelation measurement can be used for detecting deterministic 
data which may be hidden or obscured in a random background. 

As an extension of the autocorrelation function, the power 
spectral density function describes the frequency composition of the 
data in terms of the spectral density of its mean square value. 

In Chapter IV an explicit and more detailed presentation is 
included which outlines both the analytical and digital computer steps 
involved in the power spectrum analysis as used to determine the pitch 
response of the F-100C airplane to atmospheric turbulence, This 
analysis will form a mathematical extension of the autocorrelation 
function, through Fourier transformation, to yield the power spectral 


density function. 
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AIRCRAFT DESCRIPTION AND FLIGHT TEST PROCEDURES 


The basic measurement used to determine atmospheric turbulence 
in flight is the angle of attack of the aircraft as measured by flow 
vanes. The flow angles measured are corrected for the effects of 
airplane motion, Chilton (3). In this particular analysis as applied 
to the power spectrum method, a continuous time history of the vertical 
component of the truegust velocity is to be initially determined. 
Although the lateral and longitudinal components of turbulence could 
have been measured as well, the vertical component has the most signi- 
ficant effect upon the aircraft's pitch dynamics, and in the interest 
of simplification only this vertical portion of the turbulence will be 
considered. 

The vertical component of the gust velocity, from Ref. 1, is 


given by: 


Wg= Ve, -Ve.+/Qzot +4 é i 


The equation above is based upon the following assumptions as outlined 
in Ref. 1: 


1. All disturbances are small enough to allow the use 
of the angle in radians in place of the sine of the angle. 

2. Boom flexures are negligible. 

3. The effects of variation in upwash on vane indica- 

tions are negligible or allowed for by calibration. 

In conducting the flight analysis, the average values of the 
measured parameters (denoted by a bar) are included to establish the 
flight measurements as increments from the mean values for the measured 
flight interval. Therefore, the time history of the vertical compo- 


nent of the gust velocity, from Ref. 1, is: 


- 
hy = v (&, -%)-v(e-8) +/ (Gz -Q, dt +h (6-6) ©) 


A reading interval of 0.1 seconds was established as being an 
optimum value for this analysis after considering the aircraft speed at 
which the data were obtained (300 knots TAS) and the frequencies 


evaluated. 
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All turbulence data used for this analysis were obtained by 
the author during the month of June 1967 while participating in an 
industrial tour at the Flight Test Branch of NASA‘s Ames Research 
Center at Moffett Field, California. The flight test data used for 
this thesis are only a portion of the total data obtained by NASA as 
part of a larger program to study the optimum dynamic characteristics 
and response to turbulence of a variable stability F-100C airplane. 

The aircraft used for this turbulence study and its principal 
dimensions are shown in Fig. 4, from Ref. 4. It is a North American 
F-100C which is assigned to the National Aeronautics and Space Adminis- 
tration and designated as NASA aircraft number 703. The F-100 is a 
single place, low wing aircraft, powered by a J-57 series turbojet 
engine with afterburner. The wing and horizontal stabilizer are swept 
back at an angle of 45 degrees. Low speed fiight and stall character- 
istics are improved by the use of longitudinal leading edge slats on 
the standard F-100 airplane, and good longitudinal control is provided 
by an all-moving "slab" tail. The controls are completely power 
actuated, Longitudinal stick forces are provided by a simple spring 
bungee whose leverage on the stick is constant at all speed ranges 
except between a Mach number of 0.8 and 0.94 where the leverage is 
altered by a Mach number sensor. The F-100 airplane has a 400 cycle, 
115 volt A.C, electrical system and a separate 28 volt D.C. system. 

The parameters measured in flight are given by Equation (5), 
and the type and location of the flight instruments are shown in 
Table I. A continuous record of aircraft response was necessary in 
order to provide a complete and effective record of the turbulence en- 
countered, therefore an oscillograph, instead of a photopanel, was 
required. All data were recorded on a 26 channel photographically 
recording oscillograph and developed on Ektachrome Aero film for 
analysis. As can be seen in Fig. 5, the accelerometers were positioned 
as close to the aircraft's center of gravity as possible. The flow 
vane used to measure: angle of attack was located on the top of the 
vertical stabilizer as shown. A standard pitot static head was. 
utilized for airspeed measurement. It is located ahead of the aircraft 


on the tip of the F-100's standard nose boom. 
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AIRPLANE GENERAL INFORMATION T.O. 1F-100C-2-5 
Airpiane Dimensions 


AIRPLANE DIMENSIONS 2.5 FT 


NOTE Airplane principal dimensions are taken with the landing 
geor struts and tires inflated to the carrect pressures. Weights 
given ore appraximate and are nat ta be used for weight and bal- 9.7 FT 
ance computatian ar flight restrictian. 
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Statham linear strain 
g2ge accelerometer. 
Model CS-8-350 


spring restrained rate 
gyro with a Microsyn 
signal generator 


TYPE OF INS fRUMENT LOC ATION INSTRUMENT 
OF INSTRUMENT LIHITS 








Nose boom 


0-500 kts 


Top of vertical 
stabilizer e 
Station 525 














Left link 
compar tment 
Station 152 






Wheel well on 
aircraft's 
longitudinal 
centerline 
Station 315 









Left link 
compartment 
Station 13 


Longitudinal |Potentiometer with 00 | Stick linkage | 719°stick 
otick Position| cycle A.C. output on centerline deflection 


+ 28°, ape 


Horizontal 
Stabilizer 
Position 





Linear votentiometer 


pickoff using Markite 
infinite resolution vot.’ 


Altitude Statham pressure cell Nose compartment|0=-50,000 ft 
> with strain gage pickoff | Station 55 


TABLE I: 


AIRCRAFT INSTRUMENTATION: 
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Station 75 
stabilizer 
position 





NASA AERCRAFT NUMBER 703 


This particular F-100C had been modified by NASA for use as a 
variable stability airplane in order to study the effects of the 
interaction between the turbulence and the longitudinal stability 
characterdstics. A planned variatdon in C and Ch was developed 


q changed to pro- 


by NASA engineers, and these two parameters were 
vide varying degrees of pitch response. 

From basic theory, the values of oss primarily determine the 
natural frequency of the short perdod mode and are a major factor in 
determining the response of the airframe to elevator motions and gusts. 

The stability derivative, Cc. , 1s usually referred to as the 
pitch damping derivative and is 4 critical in longitudinal dynamics 
because it establishes a major portion of the damping of the short 
period mode. This damping effect comes mostly from the horizontal 
stabilizer. The appropriate oo and CF gain settings were obtained 
by means of an analog computer matching q technique to produce the 
desired combinations of the short period frequency and damping ratio. 

For this analysis, by varying Snes and Ch » 4 resultant value 
of M, equal to -14.47 and a equal to .3068 q was programmed into 
the aircraft to produce a short period frequency of .6 cycles per 
second and a damping ratio equal to 0.10. 

The simpliféed schematics in Figs. 5 and 6 show how the variable 
stability mode of this NASA F-100C is obtained. Three modes of control 
may be manually selected by the pilot: a. Basic Mode, b. "Fly By 
Wire" Mode and c. Variable Stability Mode (includes "Fly By Wire" 
mode). 

In the variable stability mode as shown in Fig. 5, a signal 
from both the angle of attack vane and the pitch rate gyro are received 
at the summing amplifier. The output from the summing amplifier is 
summed with the pilot's input and fed into a servo which controls the 
horizontal stabilizer's hydraulic actuating valve to provide the 
desired response characteristics. 

A detailed schematic of the pilot's input as a closed loop 
system is shown in Fig. 6. A potentiometer is used to detect control 
stick inputs with a sensitivity such that the basic aircraft stick-to- 
stabilizer gearing is preserved. A torque tube between the stick and 


the hydraulic valve is normally mated firmly in the Basic Aircraft Mode. 
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In the "Fly By Wire" or Variable Stability Modes, the torque tube is 
mechanically adjusted by the system to have a greater increment of 
play. This adjustment prevents any feedback from the stabilizer to 
the stick; and allows additional commands to the stabilizer in the 
Variable Stability mode. 

C and C_ are altered by the variable stability system in 


accordance with tHe following equations: 


poawe (7) 
Cring Ea Cig OQASIC PAS Cr4e| where Ay ry a 
ee 
ay )/ 
Cig eff * ia DA SIC +iGing NE a] o 


In each example, the basic values of the F-100's stability 
derivatives are modified by the second term in the equation whose 
value is altered by the servo system. 

Since the longitudinal response of the F-100 airplane to tur- 
bulence was of prime importance, all flight tests were conducted with 
no longitudinal control inputs, however the wings were kept level, 
while data were being taken, by applying minor lateral stick inputs. 
As a normal procedure, the turbulence was entered from a stabilized 
flight condition at an airspeed of 300 knots true, which is a realistic 
cruise airspeed for the F-100 airplane at the test altitudes selected 
(13,000 to 15,000 feet). Although no longitudinal control inputs 
were desired, both stabilizer position and control stick deflection 
were checked continually when analyzing the flight test data to insure 
the validity of the measured aircraft response. 

Most of the turbulence encountered was of a mountain wave or 
clear air turbulence type. The leeward side of Mount Whitney and the 
neighboring peaks in the Sierra Nevada range were used as a flight 


test area, 
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BASIC COMPUTATION PROCEDURE 


As in any statistical approach to analyzing data and predicting 
trends, the more information obtained, the more valid the study. This 
concept was tempered by the selection of that flight test data which 
were obtained with the F-100C programmed for the most realistic short 
period frequency. Accordipgly, after screening the available test data, 
NASA Flight 26, Run 8 was selected as being a valid and optimum flight 
record for the power spectrum analysis. On this flight, fhe F-100C 
was programmed for a short period frequency of .6 cycles per second, 
It is reiterated that the intent of this thesis is to demonstrate the 
application of the power spectrum analysis to determine the dynamic 
response of an aircraft to turbulence rather than to establish any new 
information on the F-100C airplane. 

The data for this flight were read every 0.1 seconds, and the 
basic parameters as outlined in Equation (5), were recorded on stan- 
dard IBM eighty column digital computer cards. Due to the large quanti- 
ties of data involved (six inflight parameters read at 920 intervals 
each) and the numerous computations required by the power spectral 
method, a digital computer is necessary to effectively process the 
data and perform the analysis efficiently. The average values of the 
measured flight parameters, denoted by a bar in Equation (6), were 
initially determined by means of a simple digital operation, and we 
itself was computed. The power spectral analysis of the measured 
flight data was performed by a digital program utilizing an IBM 
7094/7040 direct coupled system computer for the calculations and 
graphing. 

Both the autocorrelation function and the power spectral den- 
sity of the input, Me (referred to as series 2), and the output, 6 
(referred to as series 1), were determined and graphed by the computer 
program. The desired transfer function, 2. was then computed in 
both phase and amplitude as a function of gust frequency. The transfer 
function was based upon a cross correlation and cross spectral density 
function which combined the average values of 6 and w_ over a 
common frequency interval. From the computer analysis, the phase and 


amplitude of = were graphed on a cross spectral density plot 


7g 
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versus frequency. To check the accuracy of the power spectrum method, 
the theoretical solution for the phase and amplitude of & was 
obtained by solving the longitudinal dynamic equations of motion and 
the pitching moment equation about the lateral axis. The phase and 
amplitude of the transfer function obtained theoretically were also 
plotted versus frequency, and the resulting curves from both methods 


were compared. 
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DISCUSSION OF THE POWER SPECTRUM METHOD 


A detailed analysis of the steps involved in the digital com- 
puter program is presented in Appendix A; however, a proper apprecia- 
tion of the digital computation can only be realized if the mathe- 
matical theory behind the power spectrum method is outlined in detail. 

As discussed earlier, the vertical component of the atmospheric 
turbulence, oo was established as a stationary, ergodic form of 
random data. An additional stipulation is that the turbulence should 
represent a Gaussian process. For the turbulence to be Gaussian, the 
fluctuations of the turbulence must have a normal probability distri- 
bution. References 1 and 5 stipulate that if one component of the 
turbulence, w_, for instance, has a normal distribution, then normality 
of the turbulence can be assumed. Figure 7, which is reproduced from 
Ref. 1, shows the cumulative frequency distribution of “ as determined 
by NASA from a thunderstorm traverse. The plot is scaled in such a 
manner that a normal distribution of chs will result in a straight line. 
The data were obtained from a 200 second record, read at 0.1 second 
intervals. A fitted normal distribution is shown in comparison to 
the data points plotted, and it may be noted that there is a good 
correspondence between the fitted normal distribution and the traverse 
data points, except at large values of See therefore, from Ref. 1, 
the assumption that there is a Gaussian distribution of turbulence 
appears to be reasonable. It is reiterated that the properties of 
stationarity, homogeneity, and isotropy, in addition to the Gaussian 
assumptions, are implé¢dd. These properties specify, then, an invariance 
in the statistical characteristics of turbulence with respect to the 
position along the time history of the sampling interval and the 
direction through the turbulence. 

From Ref. 6, Taylor's hypothesis that the variation in gust 
intensity existing along the flight path at any given instant will 
remain substantially the same until the aircraft has traversed the 
given body of air, is utilized. The relative length of the body of 
air along the flight path, therefore is large compared with the dis- 


tance traversed in the reaction time of the aircraft. 
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Oo Measured Distribution 


Probability i 
or —L ted Normal Distributiora 
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eID, ft/ SEC 
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Vertical Component of Gust Velocity, W9 ft/sec 


Figure 7. - Probability of equaling or exceeding 
given values of vertical gust velocity 
for thunderstorm (traverse 4). 
(Reproduced from Ref. 1) 
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With the mathematical assumptions established, a detailed 
analysis of the mathematical steps behind the power spectrum analysis 
can now be presented in conjunction with appropriate sequential 
reference to the computer program which is detailed in Appendix A. 

A block diagram (Fig. 17) showing the steps involved in the 
computer analysis and the appropriate mathematical expression for 
each step is presented in Appendix A in order to provide a complete 
picture of the entire digital program involved in the power spectrum 
analysis. ; 

The magnitude of O and eo? as computed from Equation (6), 
is tabulated at 0.1 second intervals on the computer printout and 
referred to as series 1 and series 2 respectively. The values of the 
data are read from left to right in each row in chronological order. 
The autocorrelation function of r= and the autocorrelation function 
of = are computed separately by taking the product of 6 or fe at 
time t and at time t + 7, and averaging these values over a total 
observation time, T. The resulting average product will approach an 
exact. amtocorrelation function (also referred to as the autocovariance 
function) as T approaches infinity. 


In equation form the autocorrelation function is expressed as: 


PB, (2) = Leow ff xt) x(tr Wot (9) 


7-—> 02 
where K(é) represents the value of either 6 or w_ at time t, and 
X(E#r a8 is the value of 6 or w_ after a time Mire whose magnitude 
is determined by the value of 7 ‘ 

The autocorrelation functions of e and w_ are each computed 
for progressive values of 7 , that vary in 0.1 second increments, 
from 0.1 seconds to 9.99 seconds. The autocorrelation function for 
r=) is tabulated on page 7/ of Appendix A and plotted versus time 
on page 78 also. Similarly Aisa for m,, is tabulated on page 82 
and plotted versus time on page 83 in Appendix A. 

The quantity A,(>), as defined above, is a real valued function 
with a maximum value at 7 = 0, Also A, 6 2-R,(7), and &, (o) > /Ry (7 
for all 7. The mean value of x(t) can be determined from /P, (oc} ; 


i.e., it is equal to the positive square root of the autocorrelation 
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function, as the time displacement increases toward infinity (Ref. 2). 
In addition, the mean square value is equal to the autocorrelation 


function at zero time displacement, ( 7 = 0); t.e., 4 = lag? (0) 
As an example of the usefulness of the autocorrelation function, 


a simple autocorrelogram (plot of an autocorrelation function) is 


illustrated in Figure 8, 
Ges 


LN AL 
Ht VV VV 


Figure 8, - Autocorrelogram for a Sine Wave. 
(Reproduced from Ref. 2) 


Although phase information is not available with this plot, it 
can be noted that the autocorrelation function for the sine wave will 
persist over all time displacements since it represents deterministic 
data. Conversely, an autocorrelogram for narrow band random noise, 
(non-deterministic data), will diminish to zero for large time dis- 


placements, assuming a mean value of zero. This is illustrated in 


Fig. 9 


AQ) 


Figure 9, - Autocorrelogram for Random Noise 
(Reproduced from Ref. 2) 


A combined plot of the autocorrelation function of the sine 


wave superimposed with the autocorrelation function of the random 


noise will appear as in Figure 10. 
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(7) 


Oo 


Figure 10. - Combined Autocorrelograms for 
a Sine Wave and Random Noise 
(Reproduced from Ref. 2) 


The deterministic data (sine wave) can be detected from the 
random background of the wide band noise. This is the principle 
value of the autocorrelation function in that the random and non- 
random data may be separated into their respective categories through 
analysis of their combined autocorrelograms. 

The autocorrelation function of 6 » series 1, is noted to 
diminish with increasing values of time, due to the inherent longi- 
tudinal damping characteristics of the F-100 airplane. This is a 
typical response of a damped second order system to a random input. 
Refer to pages 77 and 78 of Appendix A. 

It is also to be noted that the autocorrelation function of 
w_, series 2, is also a maximum at / = 0 and decreases to a minimum 
value for increasing values of J (pages 82 and 83 , Appendix A). 
This is the expected characteristic of the autocorrelation function 
for random data. 

After establishing each autocorrelation function, [NG and 
At » the power spectral estimates of both Pn and w_ are computed 
separately and tabulated on pages 79 and 84 of Appendix A. The 
power spectral estimate is simply a Fourier transform of the auto- 
correlation function and describes the general frequency composition 
of the data in terms of the spectral density of its mean square value, 

Le (Refs. 2 and 7). The time mean square value provides a 
measure of the disturbance energy per unit time and thus has been 
referred to as power in the study of turbulence and has been termed 


the energy or power spectrum. The power spectral density function 
specifies the contribution of the harmonic components to the mean square 
value of the random data over the specified frequency range. 
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The average squared value of the random data will approach an 
exact mean square value as the observation time, T, approaches infinity 
for a frequency range between f and A f (Ref. 2). 


In equation form: 


= 
Zz ; 
[eo INA ESTAS OMG £ ak) oft (10) 
" sos 7 eo | TT. 
a 

where x(t,f, A £) is that portion of x(t) in the frequency range 
bounded by f and f + Mf. In this analysis, of course, x(t) could 
represent either © or St For small values of 4 f, the power 
spectral density may be computed to be: 


C,. (Aat = Bf feat) — 


In precirze form; 


C,Aaf= = Ln) a table 


Ars “4 Aro AT a0 (4OT; ace (12) 
wneee G,(£) represents the power spectral density function of either 
QO or Me: The power spectral density is always a real valued non- 
negative function. 
As noted earlier, the power spectral density is a Fourier trans- 
form of the autocorrelation function and for stationary data may be 


expressed as: 


a -~f27r 7 ii 
6G, 4=2/R (Ae ze A= 4/R, (cos 27f Tabhwze 


(13) 


Thus, a Ref. 2, the mean value of x(t) is given by | 
Wipe Af Eee] and the mean square value is given by Yrs = /b; (Aoaf 
The mean Square value, therefore, is equal to the area mame the plot 
of the power spectral density function versus frequency. 

From Ref. 1, typical power spectra for the vertical component 
of turbulence in clear air, cumulus clouds and thunderstorms are 
plotted in Fig. 11. The power spectral density is in units of (fe) eeene 
per radian/ft, and it is plotted versus the reduced frequency in 


radians/ft, with an additional scale of wavelength in feet, also along 


the abcissa. Since the gust disturbances are essentially space 
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Typical power spectra of vertical component of 
turbulence measured in clear air, cumulus clouds 
and thunderstorm. 
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disturbances rather than disturbances in time, a reduced frequency, 
/S2 , in radians/ft, is normally used in the power spectral analysis 
of turbulence. The relative heights of the curves in Fig. 11 indi- 
cate the variation of the intensity of the turbulence encountered 
under the different weather conditions. The similarity of the slopes 
of the power spectra is apparent, and their magnitude is equal to 
ge » where C is a constant. The spectra show a characteristic 
rapid decrease in power with increasing frequency. | 

The root mean square value of the gust velocity is a measure 
of the intensity of the turbulermce and is equal to the area under the 
portion of the spectrum evaluated. From Fig. 11, it may be noted 
that the root mean square value of the turbulence varies in magnitude 
from 4.48 ft/sec in clear air to a value of 13.38 ft/sec in a severe 
storm. The above example is presented to describe the type of infor- 
mation that is normally available from the power spectral density 
plot of a given parameter. 

The power spectral density of 8 and e are plotted on pages 
80 and 85 of Appendix A. As expected, the power spectral density 
of or (series 2) eigiueeece with increasing frequency. The power 
spectral density plot of © (series 1) also decreases with increasing 
frequency; however, a significant mode of response at the programmed 
short period frequency of .6 cycles per second is apparent. 

The general dependence of one set of data (such as 6 ) upon 


another cys) may be described in more detail using the cross correla- 


tion function. From Ref. 2, two time histories of data are illustrated 


in Fig. 12. 
x(t) 
O ~ t 
y(t) z 
2) t 


Figure 12. - Cross¢correlation measurement. 
(Reproduced from Ref, 2) 
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Like the autocorrelation function, the average product of x(t) 
at time t and y(t) at time t + J is taken over the observation time 


t. For this analysis x(t) is comparable to ra , and y(t) is comparable 
to vhs or vice versa. As T becomes infinite, the resulting average 
product approaches an exact cross correlation function. As illustra- 
ted in Ref, 2: 
r 
R yy (D2 Lams £ [xy (tr dt 
Teo 7 (14) 

ar J ) is always a real valued function which may be either posi- 
tive or negative; however, unlike the autocorrelation for random data, 
R vy CZ, may not necessarily have a maximum at J = Q or be an 
even function only. aK J ) displays symmetry about the ordinate 


when x and y are interchanged, so that st we é ig J). 


The absolute value of the cress correlation function is bounded 
by: 


1 Rey OM" SKY (o)- Ry (oO) 


Ry (MW £ BIRO) + Ry Coy 


where. RL is comparable to the autocorrelation owe SE ale of “EE i is 
comparable to the autocorrelation =e of © , and Ry is the 
cross correlation function based on 6 andw_.. The cross correlation 
function (labeled as cross covariance) for re" and w_ is tabulated on 
page 86 of Appendix A, for both a positive 7 arid a negative J , 
and plotted versus time on pages 87 and 88. 

The Fourier transform of the cross correlation is taken to 
determine the cross spectral density function which is comparable to 
the power spectral density function for a single variable. Since the 
cross correlation is not an even function, the cross spectral density 
function is a complex number composed of the co-spectral density func- 


tion (real part) and the quadrature spectral density function (imagin- 


ary part). Therefore: 
Gry (f) = Cxy Ff) - 7 Oxy (F) (15) 


43 


where Oxy is the co-spectral density fenction and Qay is the quadra- 
ture spectral density function. 

The co-spectral density function is the average product of 
x(t) and y(t) within a narrow frequency interval between f and f + 

Af divided by the frequency interval. 

The quadrature spectral density function is similar to the co- 
spectral density function, except that either x(t) or y(t) is shifted 
in time to produce a 90 degree phase shift at the frequency f. 


The co-spectral density function is expressed as: 


AF32O To 
and the quadrature spectral density function is: 


Cyy (L)= Lon bin apr /xlet at) ylt £Latlot : (16) 


Ox yY (A) = Ls Lame mar fl otar)-y(t, a rae ) Jt | (17) 


where x(t,f, 4 £) and y(t,f, & f£) are filtered portions of x(t) and 
y(t), and y°(t,£,4 f) signifies a 90 degree phase shift from y(t,f,Af). 
The co-spectral density function and the quadrature spectral density 
function for 6 and a are tabulated on page 89 of Appendix A. 

The cross spectral density function may also be expressed in 


complex polar form as: 


~7 Gy, (F) 
Cry (f= 1Gyy le pty 


(18) 


The magnitude/G,(f)/ of the cross spectral density function is 


/Cxy (Al=\V yA + Ory (O (19) 


and the phase of the cross spectral density function © yy(£), is 


Qf) 
Oxy FAY 


equal to: 


expressed as: 







—s 


ee) 7a 


The magnitude and phase of the cross spectral density function 
are tabulated on page 89 of Appendix A and plotted on pages 90 and 
91. Again a significant peak in the amplitude of the cross spectral 
density function is noted in the region of the programmed short period 
frequency of .6 cycles per second. 
The amplitude of the transfer function, between the input, 
x(t), and the output, y(t), is obtained from the amplitude of the 
cross spectzal density function, ava: and the amplitude of the 
power spectral density of the input, G(£). The amplitude of the 
transfer function, input x(t) to output y(t), is expressed as oe 
and is equal to: 
+1 yy = Cxy (21) 
Gx 
In this analysis x(t) corresponds to we and y(t) is likewise 8 ; 
so that: 


jE f= Gas) ee 


The phase of the transfer function is the same as the phase 
of the cross spectral density function as outlined above. 

The amplitude and phase of = versus frequency are tabulated 
on page 92 of Appendix A, where the ofase of the transfer function is 
expressed as portions of a circle. The appropriate translation to 
angular magnitude; in degrees, is noted on the computer program beside 


the printed figures. 


S 


Two separate curves showing the phase and amplitude of Wg 
versus frequency were plotted from the results of the power spectral 
analysis and are included as Figs. 13 and 14. It is noted that the 
peak amplitude of = , and also a phase shift of = , occurs at the 
programmed short awnies frequency of .6 cycles per second. Although 
scatter of the data points is apparent at frequency values greater 
than the short period frequency, the amplitude peak of the transfer 
function and a distinctive phase shift may be readily determined from 


the plotted curves. 
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THEORETICAL SOLUTION FOR THE PHASE AND AMPLITUDE OF THE TRANSFER 
FUNCTION 


. In order to determine the accuracy of the phase and amplitude 
of <, as computed by the power spectrum method, the appropriate 
vores al equations of motion for the F-100C airplane were solved 
on a theoretical basis. Using the programmed values of Cine and 
Ce» as well as the standard stability derivatives for the F-100C; 
the phase and amplitude of — were computed over the same frequency 
range as that used for the power spectrum analysis. 

The following basic assumptions were made in determining the 

equations of motion and solving for =. 

1. The equations of motion were written with respect to the 
Stability axes. 

2. The aircraft's elevator was considered to be locked, so 
that there were only three degrees of freedom. The equa- 
tions were written therefore for linear motion along the 
X and Z axes and for pitching moments about the Y axis. 

3. The aircraft was considered to have a linear response to 
the turbulence encountered. 

4. No lateral/longitudinal coupling occurred. 

The airspeed was stabilized at the beginning of the run 
and the engine thrust was constant throughout the run. 

6. Quasi-steady flow was assumed. 

7. Only the vertical component: of turbulence, ee contributed 
to the aircraft's normal acceleration. 

The phase and amplitude of the transfer function between e 
and 9 were initially obtained using a simplified solution in which 
the phugoid motion was considered to have a negligible effect upon the 
short period oscillations due to the turbulence. See Ref. 9. This 
simplified method is designated as method A in this section. To check 
the theoretical results, a second solution involving all of the essen- 
tial terms included in the stability equations is included. This exact 
method is designated as method B. Both data sheets for methods A and 


B are included in this analysis together with the appropriate values of 


the stability derivatives and the important terms used in each solution. 


The results of both theoretical solutions were in close agreement. 
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The phase and amplitude of - obtained by methods A and B are 
plotted in Figs. 15 and 16. It can be noted that the peak amplitude 
and the phase shift of a occurred at the programmed short period fre- 
quency of .6 cycles per second. 

From Section 13 of Ref. 9, the longitudinal equations of motion 


are developed with respect to the stability axes: 


KQXKISé A#tGE COS O, = Tu (cos 6)U t 16 pp ORPM Gos = 68) 
AX U4 XQ et KWAK Wt X gp SE +X gp FH X 5g 5b 
ZOAxIS* WG Go S177 Oy = -7y (sin €)U- Ts pepo RPM o/77 S a 
HZ UAZEG +E WH ZW +2 36 +25 3h + Z gpk 
. uO» Or BIT; t+ Palle To Span tM Ud & MoG 
one 2” a a eee” 2% (25) 


AM We My we M, SE eM SE + Mp 3h 


These equations may be further simplified since the angle ©, is 
assumed to be small so that sin O, # 0 and cos 6, = 1. In addition, 
since there is no elevator deflection, flap deflection, or speed 

brake deflection during the run, S- de and 5B all equal zero. 

Also under the initial conditions, the thrust is assumed to be con- 
stant; therefore, T. is equal to zero, and the variation of the thrust 


with RPM, T is also zero. xX _, the variation of the force along 


y 
the X axis, Bites , the variation of the force along the Z axis, due to 
the angular pitching velocity, are negligible. In addition X and a 
the variation in the forces along the X and Z axes due to the time 
rate of change of the vertical velocity, are negligible. For this 
analysis, the vertical velocity is considered to be composed of two 
components, w and wae where w is the steady state component of the 
vertical velocity and a is the fluctuation of the vertical velocity 
due to turbulence. 

Eliminating the terms as noted above, separating the vertical 


velocity into the steady state and turbulent components, and intro- 


ducing La-Placiar notation, the equations of motion become: 
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a, © mmm CL +X W + X,Wo (26) 
WwW -US8 = ZU + ZW + Z,, Wo C27) 


S76 =MYyUt Mg. 5° *M, We + MW + M SW (28) 


Re-arranging the equations and solving for the Ke terms: 


(5 -Xu)d -XyW +96 = XW (26) 
al Meche USS AINE (27) 
mead a | —(M.,S#My)W + 5(5-Mg)6= (ag W9 (28) 

Using Cramer's rule, the above determinant is modified to solve 
for so that the transfer function is equal to: 
(S-Xu) —Xy Xw 
-Z, (-2,) Zy 


My -MeS4+M,) MoM 
GE = tO ooo 


xu) —Xw Y 
—~Zu (s-z,) “3G 
~Mau - (My $+ Mw) 5(S-Mg) 


Expanding the determinant: 


oO - Kyl Zul4eS4+ My) +My (s-2)) - 2, Els-xu)(Mp 5+ Mw) 
We _ i ll 
—Mu tw] + Mwlls-XulS-2w) — Zu Xw] ae 
A 
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where 4 denotes the original determinant in the denominator of 


Equation (29a). Equation (29b) can be further simplified to: 


2 = (nMgyr Mw)s? + L[XwZu Mw + Xw Me 
W9 a 
+ 2 eee (29c) 
A 


Multiplying the numerator of the ratio of the determinants by s, an 


expression for = is obtained, so that: 
5, 


6 = (Auli My) 2 eee 
ya\ 


£6 
42 Me XG) = Pee (30a) 
Ee ammmnmnammmnemmneeamemne pesmenememeenemml 
AA 
Method A: The Simplified Solution of the Longitudinal Equations 
of Motion 


For the simplified solution, the original determinant of the 
denominator is considered to be equivalent to the product of the 
quadratic expressions for the short period and phugoid motions, in 


Laplactin: form, so that: 


6 a (Bu Mist Mw) 34 Dee Zu is + Xo Mert Be ay Kea Mi) Myo (ha PE) J 5? 


& (94 2 Foy yy Ft PHNS?t ZF sp tspS + tu"sp) (31) 


where %,,, and $%,, are the damping ratios of the phugoid and short 
period motions respectively. And &, and W/sp are equal to 277 

times the phugoid and short period motion frequencies. The contribu- 
tion of &,, is considered to be of negligible significance in the air- 
craft's dynamic response to turbulence so that the terms containing 
i/,,y Can be eliminated. Therefore, the denominator can be approxi- 


mated by: 


2 Z 


2 
Canceling s in both numerator and denominator and replacing the 


operator with jw-, the simplified expression for a becomes: 


5 a ey el) iy eed # XKyZuMiy ¢Xw Ma + Zw My Xu My) —M wp (Xa +Zw) 
ra (orsp— UW?) + 425 sp W) (32) 
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The linear velocity along the X axis is considered to be constant so 


that the perturbations Zz M and xX may be neglected, Equation (31) 


therefore reduces to: 


8 = (Zithy -Mt) + fwlZw My + My) 


W9 (Wp -WD)+Y2 7, WUspW 


The amplitude of the transfer function (method A) is: 





and the phase of the transfer function is; 


% = %O°- tan™ 2FspWsp WJ 


(Goes _ iw 2) 


2 = went, + Mw) 
wg “Ww sp ~0s?)*4 47,* ee Ww 


(33) 


(34) 


(35) 


The expressions for the stability derivatives are given by 


Ref. 9 as: 
ey Ba 
By = LEY (Chae + Co) 


SOS Eee. 
My, = $F (Come) 


V 
“Vw = —_— (Core) 


Also from Ref. 10: 


2 


sp = -La Mo - Mee 


28 sp Wsp = Loe —M¢ ~“M<x 


where: 


= = 


tf 
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(36) 


(37) 


(38) 


(39a) 


(40a) 


and therefore: 


MG = jibe <ff Ey) ie De ea (40b) 


Ag = ~Lh a ep] Mg —Lesp | 


(39b) 


The computed values of the stability derivatives noted above, 
as well as the quantities used to compute them, are listed in Table 
II. Using the values of the stability derivatives as computed in 
Table II, the phase and amplitude of — are computed for a range of 
frequencies from .05 cycles per second to 2.00 cycles per second. The 
values of the phase and amplitude oEe ‘ as computed by the simpli- 
fied solution (method A), are listed in Table III. The phase and 
amplitude of - » aS computed by the simplified solution, are plotted 


versus frequency in Figs. 15 and 16, 


Method B: The Exact Solution of the Longitudinal Equations 
of Motion. 


The initial computational steps of the simplified solution, 
method A, as well as the same given stability equations and initial 
assumptions are implied for the exact solution as well; however, for 
the exact solution the entire determinant for Ps will be retained in 
the denominator of the —- term throughout the analysis. Starting 


with Equation (30a), om is equal to: 
7 


~MllatEu)fur 





6 = (ZwM +My Ne. 
Wg 





D4 Ly Za My thay My t2y MwrXeMy (30b) 
Zs 


where A is equal to the expansion of the determinant formed by 


the longitudinal stability equations, solved for woe 


A = (S- Xu) lS -2ZwN5-Mg)=5Uo(Miys + Mw)] 


4 2L Shilo M2)4 (MS #Mw\l-Ma [SUnX,,-G(S-B,)f (438) 
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TABLE I1: FLIGHT DATA AND COMPUTED STABILITY DERIVATIVES 
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s(5—XuNS-# w(S- M9) = $y (S-Xu NM jy 5 + Mw 4 Gee (M 2 54M) 
— 5 Xw By (S-Mg) +9 Mu (S- 21) -~5U0 X, Me (43b) 
A= 5 AGG7ere Mg + Uc My)5+ LMg ky t Zw) t+ Xu Z y-Xwlu 
— 5 Aly Ku Mg dfs LUs Wy Muy Loy Ma) + Mg Win Bu -Xy Ev) 
+G(MytZyMg Ys + 9(Zudw-MuzZw) (43e) 


For simplification, the real and imaginary parts of the numera- 


tor and denominator are considered separately and designated as: 


Ry The real portion of the denominator 

I, The imaginary portion of the denominator 
Ry The real portion of the numerator 

I, The imaginary portion of the numerator 


6 


Therefore, the amplitude of the transfer function — , is 
7 


o. ee Ry & Div (44) 
ng ote TT Ee 


and the phase of the transfer function is: 


expressed as: 


it ic -—7/orr ‘IN gn ee: (45) 
Wg Rn yo 
where: 
2M, (iC . (46) 
hey = a LM ylky tin) Kyu ys —Xy Mu ~ Bw (My —XuMy,)] 
Pai” A (46a) 
ap) = Ce Lt My % My J (47) 
: ) ey) 
se aod (47a) 
Lc = Ie? fot AL iy ra Mg F PS t Us (Nu My - Xe My) 
e Hg (Hv Br Ka tw) # (Mut Bu Mn) p (48) 
7] pe Le (ae: z’) (48a) 


ao 


FD = oe wil Ce pEy)t ka Biv ead 


-~% (Mw- Ka Mw)f + Q(ZuMyw - Mua Z/) (49) 


Riz wt-w*t +O (49a) 


The stability derivatives as listed for method A, and their 
values as tabulated in Table II, apply to method B as well. In addi- 
tion, the following derivatives, computed in accordance with Ref. 9, 
are used to solve for the coefficients A through F, of Equations (46a) 
through (49a): 

Xu = oe Cr (50), Zar, ee a =e Cy (52) 


V Es 


hy = BELG-&) (5D, My =-( GE) 9 


The computed values of these stability derivatives, as well as 
the quantities used to compute them, are listed in Table II. The values 
of A, B, C, D, E, and F are listed in Table IV. The phase and amplitude 
of = are computed and listed in Table IV for a frequency range from 
.05 cycles per second to 2.00 cycles per second. The phase and ampli- 
tude of » aS computed by the exact solution, are also plotted versus 
frequency, together with the values obtained from method A, on Figs. 15 
and 16. The two different theoretical methods produced almost identi- 
cal curves. For both methods the peak amplitude of Pas » as well as the 
phase shift, occurred at the programmed short period frequency of 0.6 


cycles per second, 
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DISCUSSION OF RESULTS 


Comparing the power spectral solution and the theoretical solu- 
tion (Figs. 13, 14, 15, and 16) by matching the phase and amplitude 
curves for a , it is noted that in both cases the phase shift of the 
transfer function is in close agreement and the greatest amplitude of 
7 » as established by both solutions, occurs at the same frequency, 
0.6 cycles per second, however a discrepancy in the value of the peak 
amplitude of & is apparent. The amplitude of the transfer function, 
as determined by the power spectral method, is consistently less, by 
an almost constant difference, than the amplitude as established by the 
theoretical solution for the same range of frequencies. The validity 
of the transfer function amplitude as established theoretically is con- 
sidered to be accurate inasmuch as two different solutions produced 
identical results and the peak amplitude of 0.3 is a realistic value 
for an aircraft such as the F-100. Refer to Myslett and Smith (10). 
The power spectral solution, as derived, would imply that in flight 
the aircraft was not as disturbed by the turbulence as the theoreti- 
cal solution predicted. This is more obvious when examining the auto- 
correlation function for 8 as computed and graphed by the digital 
program. Although the distribution of the fluctuations of the auto- 
correlation function for Pa was in agreement with theory; the magni- 
tude of the fluctuations was extremely small with an average absolute 
value of .000001. 

The computational results and the observation mentioned above 
led to a close scrutiny of the instrument sensitivity on the test air- 
craft used to gather the data. [In particular, the sensitivities of 
the instruments measuring pitch angle, pitch rate, normal acceleration 
at the C.G,, and angle of attack were checked because these parameters 
provide the most vital information regarding the actual dynamic response 
of the aircraft to turbulence. The entire flight record for Run number 
8 was re-examined for the typical inflight measurements of the para- 


meters mehtioned above. The results are summarized in Table V. 
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Table V: Instrument Sensitivity and In-Flight Measurements for 
Flt. 26, Run 8, 



















Parameter Measured |Inst. Sensitivity 
(in terms of trace 
deflection) 


Average Deflection 
Flight Record 


Typical 
Resultant 
Value 






























12°/inch .22 inch ag 16 
9.1 07 inch t 7° pene 
sec/inch 
, + 
4g¢'s/inch .09 inch - .,36g 
4.5°/inch .15 inch 9° (4%) 
-0O5 inch ay: (- ey) 


The tabulations above indicate a generally small response of the 

aircraft to the clear air turbulence encountered, The data from 
Table V, the diminished magnitude of the autocorrelation function for 

6 , and the consistent trend of the power spectral solution for the 
amplitude of = to be less than the theoretical amplitude of Pe 
indicate that ae sensitivity of the instrumentation aboard the test 
aircraft was insufficient. This conclusion was further verified during 
the investigation of the instrument sensitivity at NASA when discussion 
with the NASA engineers revealed that the aircraft used in collecting 
this data had been originally instrumented for the purpose of collecting 
maneuvering type data. The instrument sensitivity requirements for this 
type of maneuvering flight are not as stringent as those necessary for 
the power spectral type analysis of turbulence. The accuracy of the 
response is even more critical in a power spectral analysis since the 
associated digital computation is of such length that any initial dis- 
crepancies will be magnified by "rounding off" errors which are carried 
throughout the computer program. The amplitude discrepancy of the 


transfer function noted above, however, was of indirect value in an 


academic sense in that it generated a thorough secondary analysis of the 
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flight test data and a better appreciation of the entire flight test 
program from considerations of the instrument location and sensitivity, 


to the digital program itself. 
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SUMMARY AND CONCLUSIONS 


Despite the single discrepancy in the solution for the transfer 
function, c , it is felt that this analysis did represent a signi- 
ficant and unique application of power spectral methods because the 
transfer function between a random input and the resultant aircraft 
dynamic response was established on a continuous basis over a finite 
frequency range. The amplitude discrepancy between the theoretical 
and power spectral solutions for the transfer function is considered 
to be entirely due to improper instrument sensitivity and not to the 
mathematical computations or the power spectral analysis itself. A 
closer correlation between the amplitudes of the transfer functions 
would be realized using the same analytical procedures if the in- 
flight response measurements were more accurate. 

Future effort in turbulence response analysis, particularly with 
power spectral applications, will require that more sensitive instru- 
mentation be incorporated on the test aircraft. In addition, the angle 
of attack vane should be re-located from its present position on the 
vertical stabilizer to a position at least two and a half fuselage 
diameters ahead of the aircraft as proposed in Refs. 1 and ll. 

Although this analysis considered a singular application of 
the power spectral method, additional areas of investigation related to 
gust loading would be worthwhile for future research. A preliminary 
study, by Press and Houbolt (12), has been conducted to determine the 
variation in normal acceleration with center of gravity position for a 
fighter type aircraft. Through power spectral analysis, the ratio of 
the root mean square values of the normal acceleration, determined by 
the area under the gust spectrum, was readily established between center 
of gravity positions. It was determined that as the center of gravity 
position is shifted aft, the peak frequency response moves to a lower 
frequency and diminishes slightly. Other significant areas of future 
study and investigation include the effect of low short period damping 
on the gust loads of aircraft in continuous turbulence and the magni- 
fication in gust load stresses that results from wing flexibility. 
Refer to Press and Tukey (13). 
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In general, the power spectral methods of analysis may be re- 
garded as an effective and accurate tool for analysis of a random type 
disturbance. In particular, it is a powerful method for further ex- 


ploration of the science of gust loading. 
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Appendix A 


The Sequence of Digital Computer Steps in the Power Spectrum Analysis 


The mathematical steps outlined in Table VI explain the compu- 


tational procedure behind each of the digital computer steps in the 


power spectrum analysis (Ref. 14). All information in Table VI refers 


to Fig. 17, and each step in the table is related to the appropriate 


computer operation by the computer sequence number. 


The symbols used in Table VI are defined as follows: 


oo value of discrete, equi-spaced time series X(t) 


ae value of time series X(t) after subtracting mean 
ee value of time series Y(t) after subtracting mean 
Lag 

Number of discrete data points 


Maximum lag 


Constant time interval 


(P) _ p (pat) 
x 


Autocovariance of series X at lag p, RY 


Power spectral estimate of series X at frequency nde 


mat 
(h) a h 7 
ES Pmat 


Smoothed power spectral estimate of series X at frequency 


h 7 
mat 





Cross-covariance between series X and Y at lag p 


Cross-power spectral estimate at frequency h /7 


mA t 
(h) 
Cospectrum, the real part of P at frequency h 77 
xy mAt 
(h) 
Quadrature spectrum, the imaginary part of Fens at frequency 
h 77 
mAt 


Smoothed cospectrum at frequency h 77 
maAt 


69 


sq6h) Smoothed quadrature spectrum at frequency h 7 
xy —— 
mat 

ane Amplitude of cross spectrum at frequency h 77 
mdt 


pHas’™? Phase of cross-spectrum at frequency h 77 
mat 


7(h) Transfer function at frequency h77 
sa mat 


TAMS Amplitude of transfer function at frequency h 77 
madat 
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TABLE VI: Sequence of Digital Computer Steps In The Power Spectrum 
Analysis 






MATHEMATICAL OPERATION REFER TO PAGE: 













The mean of the data is first 
calculated and subtracted from each 
value of X'. 


(Original Data); 










/ 
Ax = 7=p =e 


P,P = R, (pat) 


The raw estimate of the power spectrum 
is obtained by 





A) zat < NL! ad | 79 
fz Re es 
Jog — Fr 2 phn COS 177 ) A Oa. F 


where Gg O< ps7 
Fae Yo pP=O, S47? 84 


wr: Atay) (EE) 


The raw estimates of the power spectrum 
are then smoothed by "hamming" 


= ea take F 





© 


| SP, = oe 6 8 
! tw) he) | 
SR =.23R > 542+ 23P 9 chen 
7 (77 (77 -/ | 
| 7 54p 4 agp fhe 8s 


~~ 
LRT 
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TABLE VI (Continued) 






REFER TO PAGE; 






MATHEMATICAL OPERATION 





A check sum is computed to check the 
computations of the estimates. 


(0) (h) 
ons fig of Eo 


This should equal R(° the autocovari- | 
ance at zero lag. 













(6) The cross-covariance is computed by the 
formula 
f)-P 
Ry wR § Sep ; 6) 2 86,87, 
88 
Bee bs, 
xy ie e =) i jb? P=O1,2...,/77 
(A 
The crowg—spectrum, 72 ,» is given by 
h) (A) (A) 
= Cry gee 2 Oxy 
(7) | The cospectrum is obtained rom 89 
(A) vt F (A) —A 
yy = S26 + co sp" 
7 Tec plRey Rry d 
The quadrature spectrum is oalh from } 89 


Ony (A) _ sty Gly” Ry )s ”) ee 


li= Qe ee 


hk fP- O, f7) 
e es 
i (oes O< P< 


| Q) Both the cospectrum and quadrature 
spectrum are smoothed by "hamming" as 
shown in step 4, 
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TABLE VI (Continued) 


COMPUTER 


SEQUENCE 
NUMBER MATHEMATICAL OPERATION REFER TO PAGE: 





The amplitude and phase of the cross- | 89,90 
spectrum are computed. 


AM ® lsc OT [soe] 





HAS, ae ARCTAN 304,) 











5c) 89,91 
The phase angle computed is the phase 
shift of series Y with respect to 
series X. 
(12) The transfer function from series X 
to Y is given by 
(4) 
th) (A) (4) £PHAS 2 
i meen TAM », © 
x & 
xy 92 






The phase shift from X to Y is the same 
as that of the cross-spectrum. Both 
the transfer function amplitudes from 
X to Y and Y to X are printed. 
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Se er See es ee a eee lt ny pe 


SOO CLESCOHSSEOOH CEASE BE REED 


INPUT OATA TO BE PRINTED OUT---YES 















INPUT SERIES TD BE PLOTTED DUT—-- Steies L= @! 


SeR/s 2 bh 


a 
VALUE OF CONSTANT C USED IN THE PREWHITENING YRANSFORWMAYTION ZUYIeX(TeII-CXdT) --- -O 


NUMBER OF SERIES = 2 __ 






__VEDRENOING== =e 





NUMBER DF OATA PDINTS = 9320 ' 
NUMBER OF LAGS CHOSEN = 100 
NUMBER DF SELECTION CARDS = 1 





USE PREVIOUS DATA--- 





NUMBER DF VARIABLE FORMAT CARDS = 1 


() ORIGINAL OATA OF SERIES 1 


-0.00061  -0.00048 
0.00027 0.00027 
0.00052 0.00102 

-0.00023  -0.00023 
0.00140 0.0Cc077 

~0.00023 0.00115 

-0.00023 ~0.00023 
0.00127 0.00265 

-0.00023 -0.00023 

-0.00073 ~0.00086 
0.00428 0.00278 
0.00152 0.00152 

-0.00011 -0.00161 
0.00516 0200667 

-0.00249 -0-00124 

-0.0003@ -0.00036 
0.00729 0.00842 
0.00090 0.00090 
0.00064 0.00115 
0.00027 0.00027 

-9.00023 0-0COL4 

-9.00237 ~Q-003)) 

~0.00086 -0.00086 

-4.00099 -0.00099 

—9-00324 ~-0.00487 
0.00315 0.00039 

-0-00161 -0.00099 

7000262 0400249 

-~0.00086 “0.00136 

-0 .00099 -0.00099 
0.d0014 0.00039 

-0.00267 -0.00267 

~0.00199 -0.00199 
0.00366 0.00177 

-0.00161 ~0.00086 

-0.00023 ~0.00023 

“9.00149  -0.00149 

~0 .00563 -0.00487 
0.00039 0.00115 

0.00253 0.00278 

-0.00061 -0.00061 
0.00115 0.00240 

-0.00124  -0.00099 

~0.001 74 ~O0.00174 

-0.00099 -0.00099 





—=YVe e e e a ° “Ve —Ve —=VWe 
~0.00061 0.00115 0.00290 0200290 0.00278 0.00202 0.00052 ~0.00073 +*0.00324 -0.00475 
“0.00588 -0.00525 -0.00274  ~-0.C0086 0.00082 0.00077 0.00265 0.00428 0.00478 0.00441 
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~0-00136 ~0.0C06L 0.00039 0.00140 0.00315 0.00453 0.00453 0.00453 0.00453 0.00453 
0.00391 0.00165 0.00014 G.00014 - - 
0.00866 0.00491 0.00315 0.00090 -0.00011 ~0.00136 -0.00349 ~0.00668 -0.0083S -0.00625 
-0.00475 -0-00186 0.00014 c.00090 0.00240 0.00240 0.00240 0.00240 0.00202 0.00077 
_ 0200062 0.00014 0.00014 0.00014 0.00014 0.00014 0.00014 ~0.00036 -0.00011 0.00165 
0.00328 0.00453 e roe e e e 6 ~—VUe —wYe “Ve 
-0.00663 ~0.00575 ~0.00237 -0.00124 -0.00048 0.00039 0.00169 0.00240 0.002645 0.00265 
my 0.00202 9.00052 -Q.00011 -0.00086 ~0.00287 ~0.00287 ~0.00287 0.00287 ~0.00136 -0.00048 
-0.00036 0.0c027 0.00303 0.60278 a 6.00253 e e e e =Ve 
-0.00124 -0.00149 -0-00149 -C 00149 -0.00149 ~0.00073 -0.00048 -0.00023 0.00077 0.00177 
0.00190 0.00LL5 0.00027 ~0.00073 -0.00073 -0.00073 ~0.00048 ~0.00023 0.00014 0.00077 
e Oe G.00Tt5 7 e e -0.000IT ; -0.00023 —Ve “Ve =VUe 
-0.00048 ~-0.00048 ~0.00048 -0.00048 ~0.00048 0.00048 -0.000 36 -0.000 36 -0.00024 0.00027 
__ __ 0.00102 0.00127 0.00140 0.00102 0.00027 0.00036 -0.00073 ~0.00174 -0.00124 -0.00061 
0.00052 6.00177 0.00290 ~ 6.00376 — o. ° eo -O.0001T “Ve c =Ve A 
~Q 00161 ~0,.00099 -0.00023 0.00077 0.00228 0.00378 0.00879 0.00729 0.00654 0.00504 
0.00240 -0.00011 -0.00149 -0.00375 -0.00337 -0 .00299 -0.00237 ~0.00174 -0.00073 -0.00048 
~O. 6. -O0.00061 -G.001TE -0.00262 > “=0.00Z6Z <-0.00Z62Z -TU.007ZEZ “=U. =0.00Z62 
-0.00262 -0.00136 -0.00073 ~0.00036 0.00177 0.00391 0.00491 0.00491 0.00491 0.00441 
0.00240 0.00064 -0.00011 -0.00061 -0.00099 -0.00L74 -0.00174 -0 00186 -0.00186 -0.001686 
-6.00166 -0. 00186 ~O.00166 = * = ° —Ve -0.00061 —=Ve e e 
0.00014 0.00002 ~ 0.00036 ~0.00086 ~0.00086 -0.00086 -0 00099 -0.00099 -0.00099 ~0.00086 
-0.C0086 -0.00023 ~0.00011 0.00090 0.00140 0.00115 0.00039 ~0.00011 -0.000728 ~0.00086 
~=0.00149 -0.0 =VUe & ~-0.0d0186 -d.d616I -0.00073 6.00152 0.00328 C.00391 O.00391 
0.00378 -0.00011 ~0.00023 -0.C0036 ~0.00036 -0.00061 -0.00061 ~0.00036 -0 00036 0.00027 
0.00077 0.00077 Q.001L27 0.00190 0.00152 0.00090 0.00064 0.00014 -0.00011 ~0.Q00011 
-0.00011 -0.00023 -0.00036 ~Q. =VUe ° ° : e 7 e ° 
0.00014 -~0.00011 -~0.00023 -0.00036 ~0.00036 ~0.00036 -0.00048 ~0.00061 -0.00061 -0.00061 
-0.00073 -0.00086 - 0.00136 -0-001411 -9.00061 -0.00023 -0.00023 -0.00036 -0.C0086 -0.O00L11 
=0.00161  -O.00174 ~=O.0C174 -O.00LT4 —“Ve =VUe o e e H ° 
0.00391 0.00265 0.00140 -0.00036 ~0.00086 -0.00I1 74 ~0.00237 ~0.00249 ~0.00262 ~0.00262 
-0-00211 ~0.00136 -0.00061 -0.00096 0.00014 0.00064 0.00140 0.00290 0.00315 0.00290 
0.00215 -_ e 6 e & C.000l4 e “Ve “Ve - e —Ue —Ve 
-0.00048 -~0.00099 ~0.00C86 -0.00099 -0.00099 -0.000 36 0.00077 0.00240 0.0034C 0.00340 
0.00340 0.00353 0.00253 0.00215 0.00115 ~0.000L1 -0.00023 ~0-00086 ~0.00161 -0.00174 
-0.00161 ~ e “-O.00149 ~Ue Ve =Ve -0.00149 -0.00086 -0.00054 —=0.00099 
-0.00099 -0.00099 -0.000L1 0.00014 -0.00011 -0 00011 -0.00011 ~0.00011 ~0.000L1 -0.000L1 
-0.00011 -O.CCOL -0.00011 0.00014 0.00002 0.00002 0.00002 0.00002 -0.00011  ~-0.00036 
-0-00174 =0.00324 ~0.00324 -0.00237 -0.00099 0.00039 0.00102 0.00115 0.00240 0.00240 
0.00240 0.00140 C.00CL4 -0.00036 -0.00099 -0.00086 -0 .00036 0.00002 0.00014 0.00002 
-0.00023 -0.C0048 -C.00C73 -0.00061 -0.00061 -0.00061 -0.00061 -0.00061 -0.C0061 -0.00061 
-0.cO06L -0.0C06L -0.00C6L ~0.00046 -~0.00036 -0.Q000L1 -0.00011 -0.000L1 -0.C0023 -0.00073 
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0.00027 0-00027 0.00027 0.00027 . 0.00027 0.00027 _—-0.00027 _—i0.00027 
0.00002 0.00002 O« 02 200002 0.00002 0.00002 6.000028 0.00062 
0.00202 0.00290 0.00315 0 00240 0.00152 0.00064 0.000028 ~0.00023 
-0.00023 -0.00023 0.00064 0.00190 0.00315 0.00403  _—-0.0037@ ~—0.00278 
0.00014 -0.00023 -0.00L]11 ~0.00211 ~0.00211 -6.00211 -0.00149 -0.000484 
0.00115 0.00115 0.00190 0.00177 0.00090 0.00002 “0.000124 ~0.00011 
~0.00023 -~0.00089 -0.00136 -0.00249 -0.00249 -0.00249 ~0.00046 0.00064 
0.00202 0.00152 0.00077 0.00039 -0.00023 -0.00023 ~0.00029 -0.00023 
-0.00023 -0.00023 -0.00023 0.00052 6.00002 -0.00099 -0.00086 -0.00023 
“Q.00161 -0.00241 -0.00199 -0.00199 -0.00161 -0 -00136 -0.00012 0.00240 
0.00278 0.00215 0.00140 0.00039 0.00039 0.00039 0.00064 0.00115 
0.00152 0.00 340 0.00416 0.00416 0.00654 0.00654 0.00424 0.00165 
_70.00349 -0.00387 -0.00387 -~0.00362 -6.00036 0.00039 0.00177 0.00403 
0.00667 0.00604 0.00283 0.00048 “0.00211 -0.00437 -0.00437 -0.00437 
C.00014 0.00014 0.00014 0.00014 0.00014 -0.00061 -0.00061 -0.00036 
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(SECOND) OF SERIES 1 
a “=0. 0.000004 
0.-1000000 0.000004 
fee ee 0.2000000 0.000003 
0. 3000000 0.000002 
0. 4000000 0.900001 
= 0.-5000000 0.000000 
0.7000000 -0.Q000001 
0. 8000000 -0.000002 
aa = 0. 9000000 — “Ve 
1.0999999_ ae -0.000000 
1.1999999 0.000000 
1. 2999999 0.000001 
1.3999999 0.900001 
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Se Eten 429999994 = SR OOG000 
5.0999994 0.000000 
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_ 661999990 0.000001 a 
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FA POWER SPECTRAL ESTIMATES OF SERIES 1 PLOTYEO IN A LOG SCALE AGAINST FREQUENCY (CYCLES/SECOND) 
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1.550 . = 7 Ee fa re 4 
1.600 . 8 ° 1.600 
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1§.21282 13.77671 13092760 13.90212 9.42284 7.44111 0. 44024 6294663 11.28884 
12210280 7.97936 8.61612 10.23179 13.18778 15.77438 16.41386 16 .06433 15.64834 15.29233 
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10.79523 12.31425 13.63C02 13242154 13.03320 14228858 10.83886 10.36048 10.45024 10.23958 
9. 40230 6.13304 12.50450 14.9467] 13.366006 12.53961 11.44370 10.19699 1i06.i6318 ~~ 10.335135 
10.84246 10. 39247 10.04213 11216625 10.88072 8.79745 6.04909 5.08421 7.36475 9.45424 
11.79835 13.€3446 11263125 11.36029 11621494 6.27198 8.02868 9.13076 10.16107 11.01480 
11.035 34 11216100 10694672 11.-22201 11.22161 10.02038 10.20139 11.39026 10.342928 4.69062 
2.10471 7.25508 9.50407 7.43672 7.66781 8.70688 7.49752 7268425 7.61549 7255884 
9.12649 10.86250 11.75197 13246063 13.10821 10.70411 10634472 4292193 6.56989 4.52726 
12.61522 = 13.7430) = 13673772) = 16619705 = 21606220 19.47628 §=— 18015808 §=616.5377H | «6S. 14948 )~=—(15.38765 
14.63116 13.47012 11. 84286 9.89484 7.30000 3.50035 2.71851 4.47925 6.97646 4.93545 
6.96720 8.00693 7245808 §.45430 1.86626 343706 1.45160 (0666441  $496300 10694904 
11. 63239 12.35054 11.63731 11.57269 9.21297 12.84823 13.77739 15.02490 16.27242 16.16575 
13219590 12.36017 11.42819 8265964 10.79142 11292395 10.42804 9.80625 9.55196 9.16695 
9.76589 8.66221 7.50813 7064513 7.94088 9.21084 9.11178 10.15373 10093676 11.57769 
10.75781 11.59329 11.76998 1244009 12.74404 12.73400 12.72397 12.23332 11.89726 11.88963 
12.53407 13.4C637 14.52535 14.42635 14.13142 14.45963 14.52525 14617078 13.46358 13.07705 
12.80174 10.94728 11.10842 10.29630 11.74354 11.56339 11.54261 11.46760 12.50092 12.68706 
12.20654 14.25188 15.28302 8262409 7.10040 7.09245 8.05601 7.12306 7.886431 10.16580 
9.02841 12242090 12.97416 4.58978 5.80362 7.50307 7.42992 2221389 2.06170 5.77683 
4.C9788 4.10334 5.82768 §.51232 6.13178 §.98401 5.81178 § 64280 5.62821 2016622 
2.23986 2.94099 22645623 3.00813 2.64823 3.10012 4.01081 3.23404 1.95280 2.04302 
2.97298 2.29115 0.89952 0.88146 2.76986 2.74784 3.36504 4.30064 4.27862 4.25660 
7.66486 6.42442 2.52312 2282069 5.01006 6.46455 8.24171 10.31553 10.12194 9.35582 
9.88404 10.19054 7.09435 6.61671 497482 §.26411 6.76164 6.18004 6.43909 6.27486 
9.15361 7.64339 6.52691 6.51899 5 265466 3.09716 4.73765 6.21600 5.49785 7.00641 
§.31958 5.31974 5.31991 § 232971 4262370 5.09163 5.85837 5.46741 4.59657 4.29997 
5.29049 5.02686 2.81128 2240527 2.20989 2012974 2.03780 2.151C9 1.88741 1.87976 
0.74866 0.39408 -0.21870 ~-1.55230 -3.09669 -3.45793 -3.45975 -0.70336 2.57968 1.79982 
1.05706 0.42898 -2.03297 -C.02930 0.29118 -1.89379 ~1.90032 -1.89721 ~4.63511 -6.52237 
-6.11569 ~5.C7483 0.42722 -0.03672 -0.71626 0.54184 1.79754 2.72782 3.19889 3.82543 
-1.48380 -1.09507 6.21192 §.51204 4.26385 3017571 4.83964 4.68648 4.85317 5.96537 
5.92871 5291433 4225136 4.25331 4.26248 411123 4211799 6.71686 9.42935 12.05823 
10.28951 8.15274 6.C5476 3215848 4.87175 8244544 6.11253 5.71297 6.07033 7.31816 
10. 46128 9.21617 9.22940 6.73487 479002 5.09665 5.51501 5.68552 6.LLLIC 7.16098 
5.99795 §.32348 6.07614 7.56941 6.53230 5.62648 5.03656 4.89385 4.91108 4.45446 
4.62964 4.96278 5.29591 5.62664 5.64146 5.49593 7.72848 9.31490 5.33454 9.36141 
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AUTOCOVAR LANCE 

















LAG 
(SECONO) OF SERIES 2 

0. 1000000 134.240084 
0-2000000 = = 1312581932 
0.3000000 129.535526 
0. 4000000 27 w921975 
0.5000000 m —— 126-609322 
0. 6000000 125.920371 
0. 7000000 125. 466989 

- 0.8000000_ 124.593373 
0. 9000000 123.793691 
1.0000000 123.0725 32 
1.0999999 a 1S SSE ESE 
1.1999999 120.731608 
1.2999999 119.777925 

=a ele 3999999 ae NY S901 0. 
1.4999999 117.890696—_.. 
1.5999999 1176133982 
1.6999999 116.442049 
1. 7999999 115.699091 
1.86999999 114.813326 
1.9999999 114.211069 
220999999 113.906106 
221999999 113.50192T 
2.2999999 ; 13036245 
2.3999999 112.740736 
224999999 112.3342 32 
2.5999998 LLL-902419 
226999998 111.688000 
2.7999998 111.562537 
2.8999998 111.016543 
2-9999998 TIO. FEGOTS a 
320999998 109.921986 
321999998 109 .'528613 
3.2999998 7 I09.178701 a 
3.3999998 108-661991 
324999998 108.327072 
3.5999996 — | = ~~" [OT I455TF 
3.6999998 107.695281 
3.7999998 107.496702 
3.8999998 “LOT. 057416 
329999998 106.584144% 
420999997 106.180769 
4.1999997 "ao.  aLoS.sl7z2se~ _ 
422999997 105.476593 
423999996 105. 365089 
, =. = = 42559956 ° °  #« ©. : 

425999995 ' 1052246566 
426999995 104 .829396 
427999995 oa ~ - £04.419641 
428999994 104.152740 
429999994 104.060753 
500999996 F550 859989 
5.1999993 103.589747 
522999993 7 103.218030 
523999993 =a Ne : ane 
54999992 102.118905 
5 .5999992 101 .688960 
526999992 a7 a > ~~ 101.382447 a 
527999991 101. 394451 
528999991 - _ _101.581653 
5.9999990 101.319965 
6-0999990 100 .600062 
621999990 a, _ 994943790 
642999989 99.288632 
6.3999989 98585267 
6- 4999989 98001596 
65999988 972774487 
6.6999988 972353745 
6-7999988 962929824 
6.8999987 ay 96.771644 
6-9999987 962674341 
720999987 96.245048 
7.1999986 95.7a1318 
722999986 95 -475562 
723999985 952122528 
724999985 94.918305 
725999985 94.807581 
726999984 94.825804% 
727999984 94.575465 
728999984 94.218449 
729999983 93 .865534 
8.09999682 93.277043 
8.1999981 92.947791 
8.2999980 92.765900 
8.3999979 92.383802 
824999979 91.735512 
8.5999978 91172016 
8.6999977 90.622938 
8.7999976 89.914673 
8.8999975 89.238074 
8.9999974% 88.793876 
9.0999973 88. 343747 
9.1999972 87.748139 
9.2999971 87.221180 
9.3999970 86.829659 
924999969 86.542577 
Je2999960- =< ee “8529694273 
926999967 85.175899 
9.7999966 34.4601 37 
9.8999965 83.7171 33 
309999964 6323365602 
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(2) “GRAPH OF AUTOCOWMAIANCE FUNCTICN OF SERIES 2 PLOTTEQ AGAINST TIME UP TO A LAG OF 10.0000 SECONO 
































a 652542. TC GRC S25 a aa eee 
74.439 92.23% 110.028 127.823 145.617 
HroccetecceteccetessetaccetesceteccetescetecceMocectesccetecccetecscetegceeteccctescestecweMeceeteceetesest® 
-O0. + - Si, [ai Se) =e 3 -O- cat 
0.100 . * 0.100 
0-200 . ® ° 0.200 
0.300 . - : « a ___ 0300 = 
0.400 . ® ° 0.400 
0.500 + * ¢ 0.500 
0.600 . 5 ° 0.600 
0.700 . al a 
0.800 . wl ° 0.800 
0.900 . © e 0.900 
1.000 + _ ~ _ # + T. 600 
1.100 . 4 ° 1.100 
1.200 . e 1.200 
1.300 . een - 4.2300 
1.400 . ® ° 1-400 
1.500 +¢ + + 1-500 
1.600 . = — es aS eS ch 1.666 
1.700 . tf ° 1.700 
1.800 . 3 ° 1.800 
1.900 . — SRS a aT (Sn “GEE or i. 
2-000 + & + 2.000 
2-100 . a 2-100 
2-200 . -. Ce © ee s  - eeeuo. 
22-300 . z ° 2-300 
2-400 . * e 2-400 
2-500 ¢ - — 7, =e - 2.500 
2-600 . . ° 2-600 
2.700 . eB ° 2-700 
2.800 . r SS 2.560 
2.900 . 6 ° 2.900 
3.000 + ¢ + 3.000 
32100 . ee . r e ° 3.160 
3.200 . * ° 3.200 
3.300 . e ° 3.300 
3.400 . . Se | i. ini y = 3.400 
3.500 +¢ + 3.$00 
3.600 . s ° 3.600 
3.700 . - a ¢ = 3.700, 
3.800 . a ° 3.800 
3.900 . : = 2 ° 3.900 
4.000 + * + &.000 
4.100 . s ° 4.100 
_____—_—s—*#e 200 ¢ ° 4.200 
4.300 . ” ' ‘ — —— 6 i a i eee ee ee a Ce ee) ee 
4.400 . a ° 4.400 
4.500 + 7 _—_ ———— ee ee 500 
4.600 . * : 4.600 
4.700 . * ° 4.700 
4.800 . * ° 4.800 
4.900 . a A 4.900 | 
5.000 + 3 + 5.000 
§.i00 . 5 e §.100 
5.200 . Ci aa * 5200 os 
5.300 . * e 5.300 
5.400 . a - ° 5-400 
5.500 + * ee +  §.500 
5.600 . * ° 5.600 
5.700 . * ° 5.700 
5.600 . * aia Yo 3 a =) 5.800 
§.900 . & = 5.900 
6.000 +¢ * + 6.900 
6.1600 e Ct ie £ aera e e 
6.200 . s ° 6.200 
6.300 . 2 Py 6.300 
6-400 . NAF ms 7 A MSocwi = 
6.500 ¢ bg + 6-500 
6.600 . * = 2 ° 6.600 
6.700 . * : 6.700 
6.800 . + ° 6.800 
6.900 . * ane : - 6.900 
7.000 + * a ae + 7.000 
7.100 . * ° 7.100 
7.200 . * 7 o 7.200 
7-300 . * e T2300 
7-400 . bd * 7.400 
7.500 + * + 7.500 
7.600 . * ° 7.600 
7-700 . * e 7.700 
7.800 . * ° 7.800 
7.900 . ~ ° 7.900 
8.000 + bd + 8.000 
8.100 . » ° 8.100 
8.200 . * ° 8.200 
8.300 . * A 8.300 
8.4600 . Re ° 8.400 
8.500 + * + 8.500 
8.600 . * ° 8.600 
8.700 . * e 8.700 
8.800 . LS ° 8.800 
8.900 . . ° 8.900 
9.000 + * + 9.000 
9.100 . ~ = 9.100 
9.200 . * = 9.200 
See0Oee ll * 0 EEE ee Ee ooo ae js 9.300 
> ssc.  ° @## °° - » oo sep 3 ~ 9460 
9.500 + * a 9.5uU0 
9.600 . » ° 9.600 
9.700 . * C 92 7U0 
9.600 . * ° 9.400 
9.900 . ad c 9. 700 
1G.0090 + € + ie. UOS 


PCocceteccotecescteccoetesentacaeetoccetescetecccetecccotacvotessetecerseteccoetecceotecceeteccseterecteccoteccest 
65.5% 64.337 101.131 118.925 136672C 154.314 
74.439 92034 110.028 127.8273 145.617 
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TREE OTFPFERENCE I 


FREQUENCY 


a = lb Oey Sanaa 


(CYCLES/SECONO) 


sls 


~ 202000 


-0. 
0.050000 
0.100000 
0.150600 
0.200000 
0.250000 
0.380000 
0.350000 


0.400000 


4 
0.500000 
0.550000 
0.600000 
0.650000 
0.700000 
0.750000 
0.800000 
0.850000 
0.900060 
0.950000 
1.000000 
1.050000 
1.400000 
1.150000 
1.200000 
12250000 
1.300000 
1.350000 
1.400000 
1.450000 
1.500000 
1.550000 
1.600000 

1.650000 
1.700000 
1.750000 


1.800060 


1.850000 
1.900000 
1.950606 
2.000000 
2.050000 
2.100000 
2-150000 
22200000 


2.300000 
2.350000 
2-400000 
2.450000 
2.500000 
2.600000 
2.650000 
2.700000 
2.750000 
2.800000 
2.850000 
22900000 
2.950000 
3.000000 
3.050000 
3.100000 
30149999 
32199999 
3.249999 
32299999 
3.349999 
3.399999 
32449999 
3.499999 
32549999 
32599999 
3.649999 
3.699999 
3.749999 
3.799999 
32849999 
3.899999 
3.949999 
3.999999 
4.049999 
4.099999 
42149999 
4.199999 
42249999 
42299999 
4.349999 
4.399999 
46449999 
4.499999 
42549999 
4599999 
4 6649999 
4.699999 
4.749999 
4.799998 
“12 649998 
4 -e-899996 
#2949998 
42999996 
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POWER SPECTRAL ESTIMATES 


O.F30079G 


OF SERLES 2 


38S. 7ZUSGET 
£82.6911449 
17.6962090 

~  Oeh&62Z8992 


5 9579135 


420517607 
Oe SET ZIT 


263770285 
1.6715009 


027405255 
0.6149712 


—_— —— —_ + 


a a a a 





0.8312680 
0-7522334 


ee 


0.8241192 
0.7903058 


0.5986931 
0. 7000603 
0.6201725 
0.5461514 
0.6196621 








0. 4137453 
023527345 
0.4807013 








0.4302773 


062442606 
0.2064999 
0.2690425 
022498153 
0.1055194 


0.1510822 — 


0-216996% 
022253814 
022266283 
0.1850663 
0.1082196 


oS. —- -— CBO. — a 


to 


SPECTRAL 


=e -uNUL6E9b 


! STEPATES JS 
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0.1065000 


021020795 
s~ 


0-1398452 
0-1294576 


0.0553440 
0.0858689 


~~ 001243207 


0.0852670 
0.0582651 


00524829 
0-0591246 
0.0817738 
0.1007028 
0-1209391 
0.1759721 
0.1336417 
0.0824%582 
0.0847758 
0.06304%63 
0.0583861 
0.0793020 
0.069926] 
00520203 
020446261 
0.0475723 
0.0300561 
0.0187037 
0.0194987 
0.0272510 
0.0330077 
0.0339645 
0.0363169 
0-0403523 
0.0334528 
0.0377382 
020427226 
0-0421150 
0-0418114 
0.0300238 
020245524 
0.0303787 
00441823 
0.0509428 
0.0358283 
0.0291407 
0.0361181 
0.0316606 
0.0365714 
0.0405647 
0.0361531 
0.0318768 
0.0233634 
0.0236066 


066334033" 


0.0361 747 
0.0302330 
0.0338973 
0.-0402399 


a 0-0681205.—~- 


136. 7196331 ANC SHOULL LE 


Nereis Hees 


ae #) POWER SPECTYRAL ESYIMATES OF SERVES 2  BLOTVYED IW A LOG STALE AGAINSY FREQUENCY (CYCLES/SECONDI _ 
-3.979 1.993 -0.007 1.978 3.964 5.950 
-2. 986 -1.000 0.985 2-971 4.95 
Hoc cetecseteccateccs steccatasseteocrectecscotercsetecceteccetdcsetaccets cootaccctessetosbetacscstesestesert 
-0. + dhs -0. ——— 
0.050 . 7 ic “a ae * ° 0.056 
0.100 e 6 ° 0.100 
0.150 . © x ° 0.150 
0-200 . — a ° 0.200 
0.250 ¢ v + 0.250 
0.300 . ¢ ° 0.300 
0.350 . = aa : a : e 5.350 ‘72 
0.400 . td ° 0.400 
0.450 . * ° 0.450 
0.500 + a _ ¢ 0-500 ~— 
0.550 . ° 0.550 
0.600 . ‘ 0.600 
0.650 . eT a . 0.650 
0.700 . e 0.700 
0.750 + + 0.750 
0.800 © 00 . * ° 0.800 
0.850 . * e 0.850 
0.900 . e 0.900 
0.950 . _ re ae ° G.950 
1.000 ¢ ¢ 1.000 
1.050 . e 1.050 
1.100 . a : ° 1.160 
1-150 . ° 1.150 
1.200 . * e 1.200 
1.250 ¢ + ae + 52505 2 = 
1.300 . * e 1.300 
1.350 . 6 ° 1.350 
1.4060 . —_ ° 1.400 — 
1.450 . * ° 1.450 
1.500 ¢ * ¢ 1.500 
1.550 . e 1.550 
1.600 . e 1.600 
1.650 . * ; : 1.650 
1.700 . - “\ . 1.700 
1.750 + + 1.750 
1.800 . _— ; —an - 1.800 
1.850 . * ; . 1.850 
1.900 . * . 1.900 
1.950 .« a oe : 1.950 
2-000 +¢ * + 2.000 
2.050 . ° 2.050 
2.100 . PON... Tae 7 <.. : 2.100 
ian mouse © eR rn a em or tr ee he e 
2-200 . 7 e 2-200 
2.250 ¢ _ a. ee + 2.2500 | 
2-300 . * 2.300 
2-350 . 2.350 
2.400 . * a. 2.400 
2-450 . 7 e 2-450 
2.500 ¢ > + 2-500 
2.550 . “ ; a = 2.550 
2-600 . * e 2-600 
2-650 . * e 2-650 
22700 « oe _— Re 2.700 
2-750 ¢ + 2.750 
2.800 . * e 2-800 
2.850 . . * _. 2.850 
2.900 . * ° 2-900 
2.950 . e 2.950 
3.000 «¢ = a + 3.000 
3.050 . C * 3.050 
3-100 . * ° 3.100 
30150 « : nn ee 
3.200 . * ° 3.200 
3.250 + * + 3.250 
3.300 e * e 3,300 
3.350 . a e 3.350 
3.400 . * e 3.400 
34450 .* a ° 3.450 
3.500 +* + 3.500 
32550 . * e 3.550 
3.600 . * e 3.600 
3.650 . * ° 3.650 
3-700 . * e 3-700 
3.750 ¢ * : + 3.750 
3.800 . * ° 3.800 
3-850 . * ° 3-850 
3.900 . * ° 3.900 
3.950 . Ld . 3.950 
4.000 +¢ * + 4.000 
4.050 . * e 4.050 
4.100 . Ld ° 4.100 
4-150 . * ° 4.150 
4.200 . * e 4.200 
4.250 + * + 4.250 
4.300 . * ° 4.300 
4.350 . * e 4.350 
4-400 e * = 4.400 
4.450 . * e 42450 
4.500 * * + 4.500 
4.550 . * e 4.550 
4-600 . * e 4.600 
4.650 . * 2 : . A 4.650 
<0. * @@=@=&« °° °° °° = 7 7 « ree 4.706 
4.750 + * + 4.750 
4.800 . » e 4.600 
4.850 . * ° 4.850 
4.900 . . ° 4.9U0 
4.950 . - ° 4.7930 
e00O - . + 2 OM 
to cce Ta csetoceatececteacodtoesseeteccetecesteccetecectececeteceet 6460s * Geiss t belts * 6 a ale ® alelalett clalcie* clchelecticielsielt 
79 UG 1.993 -0.00/ 1.976 3.904 5.95u 
—2e'1Bt -1.000 0.985 2.971 4.957 



























































. CROSSCOVARIANCE = = = 


= LAG © CROSSCOVARIANCE 
OF SERIES 1 ANO 2 


OF SERIES 1 ANO 2 





























(SECONO) (POSITIVE Tau) (NEGATIVE TAU) » 
- ~ 0.2600 4 7 e “Ve 
0.3000 -0.000740 -0.000942 
0.4000 -0.001713 -90.000726 _ 
0.5660 -02002795 -6.660261 ae 
0.6000 -0.003751 0.000370 
0. 7000 -0-004492 0.001021 
° “=6.004957 ° 
1.0000 -0.004926 _ 0.002243 
1.1000 -0.004590 e 2 6 
1.2000 -0.004196 0.002163 
1.4000 -0.003308 02001660 
1.5000 -0.002807 0.001329 
1.6000 0.002360 6.001004 
1.7000 -0.002002 0.000735 
1.8000 -0.001 768 0000464 
1.9000 -0.001724 0.000197 
2.0000 -0.001740 0.000116 
2.21000 -0.001862 0.000132 
2.2000 -0. 002088 0.000140 
Ze 5000 = -0. 002610 0.000549 
2.6000 -0.002579 0.000745 
2.7000 -0.002411 0.000938 
2.8000 ~0.002174 02001039 a 
3.0000 -0.001673 0.001123 
3.1000 -0-001521 0.001152 _ - 
3.2000 -0.001494 0.001056 
3.3000 -0.001567 0.000874 
3.4000 -0.001759 _ 02000754 Ss a 
3.5000 -0.001926 0.000637 
3.6000 -0.002033 0.000584 
_ 3.7000 —_ -0.002031 0.000659 
3.8000 -0.001997 0.000765 
3.9000 -0-001960 02000875 
4.0000 79.001925 0.001066 
4.1000 -0.001892 0.001312 
4.2000 -0.001829 0.001472 
a 42.3000. _ -0.001749. —— OT ll 
7 =a 424000 = ————si ON OT 0.601692 — ee 
4.5000 ~0.001462 0.001714 
426000 00001372 a re 1+) 
4.7000 ~0.001333 02001405 
429900 -Q.001378 ae CUR ae 
5.0000 -0.001471 0.000837 
5.1000 -0.001531 0.000829 
522000 -0.001572 gee 0 z ee 
5.3000 -0.001691 0.001045 
§.5000 -0.002099 fee «8 a = 
5.6000 -0.-002260 0.001686 
5- 7000 -0.002391 0.001966 
5.9000 ~0.002363 0.002483 
6.0000 -0.00229% 0.002483 
6-1000 -0.002250 _ ~ 92002316 rT ‘s anne 
6.3000 -0.002155 0.001937 
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GRAPH OF AMPLITUOE OF CROSS-SPECTRUM OF SERIES 1 ANO 2 (IN LOG SCALE) AGAINST FREQUENCY (CYCLES/SECONO) 
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